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Abstract. - We develop a new theoretical framework for describing steady-state quantum 
transport phenomena, based on the general maximum-entropy principle of non-equilibrium 
statistical mechanics. The general form of the many-body density matrix is derived, which 
contains the invariant part of the current operator that guarantees the non-equilibrium and 
steady-state character of the ensemble. Several examples of the theory are given, demonstrat- 
ing the relationship of the present treatment to the widely-used scattering-states occupation 
schemes at the level of the self-consistent single-particle approximation. 



During the last few years many ab-initio calculations have addressed the electronic struc- 
ture of systems with nonzero electrical current |ll2l."i| . We will refer to these as the occupation 
scheme approaches (OS) since in practise one occupies the right- and left-going scattering 
states up to two different electrochemical potentials, fj, R and /i L respectively. The essential 
idea of the OS come from Landauer's treatment of coherent transport in terms of the trans- 
mission matrix of the conductor 0]. Later, Caroli et al. |S] and independently, Feuchtwang 6 
developed a formal theory of tunnelling based on the technique of Keldysh non-equilibrium 
Green's functions [7j which could be extended to address coherent transport as well. Recently 
it has been shown |Hj that these two approaches are indeed equivalent at the level of a single- 
particle approximation in the spirit of Kohn-Sham density-functional theory. However, the 
latter approximation is very hard to justify. It is not clear what effective potential one should 
use; the use of the local density approximation (LDA) is a mere hope rather than a secure 
approximation. We also believe that a certain difficulty might lie in the formal search for 
steady-state non-equilibrium Green's functions using a unitary (Hamiltonian-driven) evolu- 
tion for t — ► oo from an undisturbed system. For example, when we adiabatically turn on 
an external field, the Keldysh technique predicts no change in temperature, in contradiction 
with statistical thermodynamics. We therefore believe that any alternative point of view is of 
great utility here. 

We build such an alternative theory using the generalised maximum-entropy principle as 
established by Jaynes [9]. Similar ideas were heavily exploited in the development of the 



(*) E-mail: pb209york.ac.uk 



2 



PREPRINT 



projector techniques for non-equilibrium statistical mechanics by Mori [lUHllj . yet detailed 
application to concrete problems is not widespread. Of the few papers, let us mention those 
of Ng |12j . and Heinonen and Johnson 13 , who consider current-carrying ensembles, similar 
to ours. However, in these papers, the essential steady-state character is not considered. This 
results in differences between the predictions of OS and maximum-entropy methods, in the 
linear-response regime, that are not present in our work. We give a coherent formalism that 
does not depend on the complexity of the system, i.e. without restriction to non-interacting 
particles or simple band-structure models. 

The statistical density matrix (DM), which represents an ensemble with known or con- 
trolled averages of given operators (Ai) — Tr[pAi], is obtained by maximising the information 
entropy S[p] — — Tr [/31og(/3)], subject to constraints on the traces of the above-mentioned 
operators 9 . In the case of quantum transport, experiments suggest that for a given tem- 
perature, composition and total current we obtain a well-defined thermodynamic state (or, in 
the case of N-shaped I- V curves, a small number of states differing by applied bias voltage) . 

Firstly, the total energy is conserved. This constraint is associated with the Lagrange 
multiplier /3, corresponding to inverse temperature for equilibrium - or near-equilibrium - 
systems. Similarly, the number of electrons is conserved and on average is given by the total 
positive charge in the background, i.e. atomic nuclei. Therefore, a constraint on the number 
of particles is used, with the usual symbol p, for the related Lagrange multiplier. The total 
current I should be the next thermodynamical parameter of the theory. On the contrary, the 
vast majority of present approaches to quantum transport use the applied bias AV^ instead. 
However, AV is not convenient for it is defined uniquely only between two ideal reservoirs, 
each being in equilibrium. These should never be a part of a practical calculation, to say 
nothing of the strongly non-local character of this quantity. In contrast, the current flowing 
through the system is represented by a simple operator and is well-defined even in the strongly 
non-equilibrium regime. We use the symbol A for the Lagrange multiplier accompanying the 
current constraint, and we later show that A is universally related to AV. 

Finally we impose the steady-state condition p, H = 0. For this to have a nontrivial 
solution, the system must be infinite along the direction of the current. Otherwise, the 
only steady state would correspond to zero current. This is equally present in the Keldysh 
formalism, where one has to consider the limit of infinite size first, and only afterwards can 
the time-evolution of the response to the turned-on transfer Hamiltonian go to infinity. To 
implement the steady-state constraint wc write the steady-state condition in any complete 
set of states (E, a\ p, H \E' , a') = for all E, E', a, a' . This particular notation stresses the 

fact that we work with a continuum of eigenstates of H , normalised to a delta function of 
energy ( 1 ). The index a runs over the discrete set of degenerate states at energy E. Each 
of these equations must be now guaranteed, with a separate Lagrange multiplier X a a i(E, E') 
and the expression in the functional to be maximised can be manipulated into 



dEdE' V ,a (E' ,E){E,a\ 



p,H 



\E'a' 



Tr 



A [p,H] 



Tr 



pL 



(1) 



where we have introduced L = [H, A]. This form is suitable for the variation with respect to 
the DM. 

Collecting all the above terms we obtain the variational condition 



S {-(log(p)) + (fl + - 0{H) + Pp{N) + (3A(I) -p{L}} 



0. 



(2) 



(^The final results obtained do not depend on the particular choice of normalisation. 
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The term (Q + 1)(1) guarantees the normalisation of the DM. We also note that we have 
deliberately introduced the parameter (3 in the definition of all the other multipliers so that 
the limit f3 — > oo can be conveniently studied. As a result of variation we obtain the stationary 
non-equilibrium DM p = exp{f2 — (3K}, where K = H — pN — AI + L. The practicality of this 
expression relies on the knowledge of the L operator. We obtain its form from the solution 



of 



p,H 



0, as an equation for L. Expanding the DM in terms of K, we see that this is 



equivalent to 



-AI + L,H 



0. If we cast the last expression in the representation of the 



eigenstates of H , it is seen that the role of L is to remove the off-diagonal elements of the 
current operator. We shall show below that L should be of the form 

La, a' (E, E') = I a%af (E, E') (a-AS{E-ET^, (3) 

where A is a finite constant, related to A as A — ire A, with e ~ \jl an infinitely small energy 
inversely proportional to the length of the system. The result is a finite number, since by 
inspection of Eq. (3) we deduce that A ~ I due to I ~ 1. 

To prove ® we set Y = —A I,H and write in the basis of {\E, a)} the equation for L 

as: 

iY a , a ' (E, E') + iL a , a i (E, E')(E' - E) = 0, (4) 

for E ^ E' we have 

L a ,a> (E, E') = iY *^v? =AI a , a ,{E,E>). (5) 

Because Y(E, E') is a result of a commutator, it is also proportional to E — E' . However, we 
need L a a *(E, E') to be zero for E = E' to keep the current at its given value and to satisfy 
the fact that L is also result of a commutator with H. This is uniquely achieved by 

L a . a , (E, E') = I a , a , (E, E') \A - A lim enS e (E - E'fj , (6) 

where 

which manifestly satisfies both conditions, since for E = E' , TreS e = 1. This is the stated 
result Eq. © can be now written in basis-independent form as 

1 <- T 



l = a[i-i»), i^y^-jjm, (8) 

where the operator I has the form of the invariant part of the current operator with respect 
to the time evolution, introduced by R. Kubo in the linear response theory where the 
timc-dcpcndcnce of the operator I(t) is determined by the Hamiltonian H. If we insert for 
the L in the stationary DM obtained from Eq. (0) the solution Eq. JS}, we obtain the final 
result for the statistical density matrix 

p = cxp{0 - (3(H - fj,N - AI Q )}. (9) 
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Fig. 1 - Compound system 1+2 (see text) and related quantities. 



This is the general form, valid even for a fully interacting system. It is an interesting obser- 
vation, that a sole requirement of the constraint on the time average of current operator is 
equivalent to the strong stationarity of the DM and the constraint on the current operator. 

We will now deduce the meaning of the parameter A. Consider we have two steady-state 
systems 1 and 2, with nonzero currents (see Fig^l described by their respective DM's Px<2)- 
The current is maintained with the parameters Au 2 ) , but we can equally well imagine, that 
far left and far right there are ideal reservoirs to which we apply bias AVi(2), so that A and 
AV is in correspondence. If we were to describe a single compound system, comprising weakly 
coupled systems 1 and 2, with only the total current being known, the DM would have had the 
form p 1+2 = cxp {(3 [Hi + H 2 - fi{Ni + N 2 ) - A(I^ + 1% )) }. On the other hand, if we weakly 
couple the originally disconnected 1 and 2, we have p\,i — Pif>2- Clearly, the averages with 
respect to P1+2 and pi^ will be the same if /3i = 02, pi = p.2 and A\ ~ A2, i.e. no change of 
total current, as well as of I\ and I2, is introduced by coupling. Exactly this happens when the 
applied biases AVi( 2 ) are identical, so we conclude that AV should be a universal function of 
A, p and f3. The universality comes from the fact that systems 1 and 2 are arbitrary. The 'Am- 
meter' could be represented by a Landauer's concept of infinitely large reservoirs adiabatically 
connected through ID conducting channel ^H- We leave detailed analysis of this situation for 
a future paper and here infer the A — AV relation from specific results in the following, giving 
AV = 2A for small A. This is a very general thermodynamical statement and removes the 
detailed considerations of near-equilibrium in reservoirs from the actual transport problem of 
interest in the nano-contact. 

We will now demonstrate several features of the general theory developed above, at the 
level of a self-consistent single- particle approximation. In the single-particle approximation, 
is sufficient to know the single-particle density matrix for evaluation of any quantity, in our 
case the current and the electron density. These are well-defined for infinite system, unlike 
the total energy or total number of particles. Since we deal with a system that is genuinely 
infinite, i.e. there is a potential drop when comparing the right and left asymptotic regions, 
with uniform current flowing, we need to resort to Matsubara Green's function techniques to 
obtain the density matrix unambiguously. The result is 



V ' ; J ^ e 0{E-AI a {E)-n) + V ; 

a Fermi-like distribution with the effective dispersion E(k) = E — AI a (E). The XE,a{x) 
diagonalise the effective Hamiltonian K. I a {E) = ±\t(E)\^/ r ^ are the eigenvalues of the 
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invariant current operator, which in the basis of right- and left- going energy normalised 
scattering states has the form 

2ttI°(E) = 

The states XE,a( x ) are unitary transformation of the scattering states given by the eigenvec- 
tors of the matrix I°(E) at each energy level E. t,r and t,f are the usual forward and 
backward transmission and reflection coefficients respectively, and finally k — \J2E and 
k = \]2{E + A</>) with A<j) being the drop in electrostatic potential energy. Crucially, the 
scattering states appear here just as a convenient complete set of eigenstates of the Hamil- 
tonian and it is the states XE.a which are actually being occupied according to Fermi-like 
occupancies in Eq. (JTDJ . In the limiting case of |r(.E)| — > we obtain the original right- and 
left- going scattering states, in agreement with the occupation scheme. On the other hand, 
for \t(E)\ — > we get nearly their symmetric and antisymmetric combinations. We discuss 
the physical significance of these in later paragraphs. 

Next we give our motivation for the identification of A as the applied bias through Ay = 
2A. We look at the expectation value of the current operator in a ID perfect wire. In the 
small A limit we have 

J=2V [ X - dE i 2 * (-ir\t(E)u[^=^2A\t(E F )\ 2 . (12) 

Since it is an experimentally well-established fact that the conductivity of a ID channel is 
2e 2 /h |15| . we can directly identify 2 A with the bias applied between two equilibrium reser- 
voirs. Due to the general arguments above we know that this relation is universal (for small 
A), so it needs to have the same form for any system. Eq. l|12f) is in complete agreement with 
Landauer's formula 0] even though it comes from rather different considerations. 

In the following we will be concerned with the self-consistent determination of the drop 
in electrostatic potential A</>, and a detailed discussion of the difference between Acj> and 
the applied bias AV. Specifically, let us suppose that our system consists of two identical 
_D-dimensional jellium-like leads. Local neutrality requires ((3 — > oo): 

fdE r / x f° dk D 

/ tt [n{x -> -oo) - n(x -► oo) = / t^Td- (13) 
Ja 2vr J_ A<p {2tx) u 

The meaning of this is that the charge appearing below the potential drop, on the right, must 
be exactly compensated by the charge transfered to the left by means of the occupancies in 
Eq. I|ll)[l . We can analytically evaluate the left-hand side for small A, obtaining 

A0 = 2A\t\ 2 \r\ 2 = AV\t\ 2 \r\ 2 , (14) 

independent of dimensionality D. Through this we can relate the 4-point conductance Gap = 
I/A<f> to the 2-point conductance G — I/AV. We immediately see, that the former gives a 
surprising result G±p = ^--pp, approaching the quantum of conductance for \t\ — > 0. This 
contra-intuitive result can be understood in terms of the occupation of nearly anti-symmetric 
admixtures of right- and left- going scattering states, present in xe.+(x). While this comes 
out of our formalism, we can expect that these combinations for the weakly connected system 
will be destroyed by a finite lifetime of the single-particle states, arising from whatever weak 
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Fig. 2 - The effective (E(k)) and the true energy (E(k)) dispersion relations with the corresponding 
Fermi energies n, jJ,+ and 



scattering by phonons or other electrons. If we model this fact by cancelling the off-diagonal 
terms in the invariant current matrix Eq. the resulting 2-point conductance turns out 

to be G = ^f-|£| 4 and the 4-point conductance G±p — ^-|£U, while the relation in Eq. I|14|) 

remains unchanged. G^p obtained without the off-diagonal terms is in complete agreement 
with the seminal work of Buttiker et al. 0] , while the 2-point formula gives the conductance 
smaller by factor \t\ 2 . We would also like to stress, that the off-diagonal elements could 
possibly play role for situations, when \t\ ~ 1 and therefore lead to higher conductances than 
those obtained from the Landauer 2-point formula. 

In order to elaborate the relation between the maximum entropy theory and OS we notice, 
that even though we work with only one parameter related to the number of particles, /z, 
from (|10f> we see, that we can define two auxiliary Fermi energies /j,± up to which the states 
a = ± are occupied from fj,± — AI± ) = fj, (see Fig. [5J • In the linear response we get 
A/Li = jit_|_ — /!_ = 2A\t\ which together with Eq. ^] results in A<fi — Afi\t\\r\ 2 . Similarly, 
without the off-diagonal elements we have A/i = 2yl|i| 2 and A<fi = A/i|r| 2 . The latter relations 
demonstrate most clearly the difference between the maximum entropy and OS. Firstly, when 
ignoring the off-diagonals, the applied bias AV in the OS is heuristically identified with 
A/i while in our treatment the thermodynamical arguments given in the first part of this 
letter suggest 2A = A[i/\t\ 2 . Second, the right- and left- going states are not unexceptionally 
occupied by the left and right reservoir respectively; we believe that particularly for \t\ ~ 1 can 
this effect be experimentally verified based on the differences between conductances coming 
from these two approaches. 

In conclusion, we have shown how the maximum-entropy formalism can be applied for 
non-equilibrium steady states. We have derived the statistical density matrix introducing the 
L-operator that guarantees the steady-state character of the statistical ensemble and identified 
its resolution with Kubo's invariant part of the current operator. A Lagrange multiplier 2A, 
conjugate to the current operator, represents the applied bias. In the second part of the paper 
we have demonstrated the theory on simple examples, discussing in detail the character of 
the density matrix within the single-particle approximation. We have shown that for systems 
with no reflection probability our theory gives results identical to the usual occupation scheme. 
This agreement slowly breaks down as the transmission is decreased, as the relevant states 
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became a coherent combinations of right- and left- going states. We have derived a simple 
dimensionality-dependent formula for the electrostatic potential drop and discussed its relation 
to the applied bias within the context of our theory. 
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